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Abstract 

^ ' The one-loop free energy of the four- dimensional compact QED, which is known to be equiv- 
^ • alent to the vector Sine-Gordon model, is calculated in the strong coupling regime. In the case, 
p^. when the norm of the strength tensor of the saddle-point value of the corresponding Sine-Gordon 
D model is much larger than the typical inverse area of a loop in the gas of the monopole rings, the 
obtained free energy decays exponentially vs. this norm. In the opposite case, when the dominant 
configuration of the Sine-Gordon model is identically zero, the resulting free energy decays with 
^ ' the growth of loops as an exponent of the inverse square of their typical area. 
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It is known [1,2], that the so-called strong coupling regime of the 4D compact QED, when its 
coupling constant g is larger than some critical value, corresponds to the confining phase of this 
theory. In the present Letter, we shall calculate the one-loop free energy of the 4D compact QED 
in such a case, i.e. in the case when g ^ 1. The expression for the partition function of this model 
has been found in Ref. [1] and has the following form 
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j Vx^{t) exp (-S^^S^^G^^) cos (s^^(9[^Xz. 
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Here ^ is an UV momentum cutoff, G^ib = where S^^-* is a dimensionless constant sym- 

metric tensor, and Sq,^ = § Xadxp is the so-called tensor area [3] associated with the contour 

c 

C : Xa(t), < t < T, Xa{0) = XaiT). Notice that T^aa'^ficrGaii is just the self-action of a single 
monopole loop, whereas Eq. (1) is the partition function of a gas of such loops, where in the 
confining phase under consideration the large loops are dominant, i.e. this gas is not dilute. 

Without loss of generality, the symmetric matrix G^p could be chosen to be a diagonal one. 
Then, since 

T^aa^p^Gap = <j> dx„ j> dy^XaVpGap = X! '^"a / ^^<^ / '^y^ ~ y^^f ' 
c c c c 
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one can see that (Gq-q) * is a typical size of a monopole loop in the a-direction (of course, all the 
four components of the matrix Gq,/3 are positive by definition). In what follows, we shall consider 
only the loops whose sizes along at least three axes, x, y, and are equal to each other, so that 
the matrix Gap has the form Gap = diag (Ai, Ai, Ai, A2). Then one can linearize the term bilocal in 
the tensor area in Eq. (1) by introducing the integration over an auxiliary antisymmetric matrixQ 
as follows 

1 r 



2V27T^Xf (Ai + A 



where 



and 
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J dB = J dBi2dBi^dBi4dB23dB24dB^4, 



Aap = — diag I 1, 1, 1, 



4Ai ° V Ai + ^2 y 

Notice, that since both Ai and A2 are positive, there could not be any singularities neither in the 
matrix Aap, nor in the denominator on the R.H.S. of Eq. (2). 
In the regime we are interested with, when g ^ 1, the term 



^Such a linearization for the case when Gaf3 ~ <5a/3 has been first proposed in [4]. 
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in the functional integral on the R.H.S. of Eq. (1) is dominant. Therefore, up to an inessential 
constant, the one-loop free energy following from Eq. (1) will be equal to the value of the ex- 
pression, standing in the second line of this equation, which it acquires at the saddle-point of the 
integral 

This saddle-point value obviously reads X/j. — ^^uF^ij,, where F^i^ is a constant antisymmetric 
matrix, which in particular could be identically zero. Hence, 



In order to calculate the integral over contours standing on the R.H.S. of Eq. (3), let us regularize 

it by introducing a small kinetic term c " , after which this integral becomes nothing else, 
but the one- loop unsubtracted Euler-Heisenberg-Schiwinger Lagrangian in scalar QED [5], and 
Eq. (3) yields [6,7] 
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sin {2g'^aiT) sinh (2g%iT) 
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where 
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i — 1,2. Here Ei,Bi and -^2,-82 stand for the "electric" and "magnetic fields" corresponding to 
the "field strength tensors" Fa/3 ~ B^p and Fap + -Bq/j respectively. Therefore we have replaced 
the integral over contours standing in the initial expression (1) by the averaging over an auxiliary 
matrix 5^/3 with a certain Gaussian weight, which is obviously much more simple. 

Let us now proceed with performing this average in the simplest case when the loops are flat, 
i.e. only two components of the tensor area (which for concreteness will be chosen to be the 
components E23 and E32 respectively) are nonzero. In this case, Eq. (4) yields [7] 
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where F = F23 is positive by definition. Since g ^ 1, in what follows we shall always assume that 

The integral on the R.H.S. of Eq. (5) can be most easily carried out in the two following cases. 
First of them is the case of the strong field x^, when F ^ -\/Ai. In this case, one finds 

g^^frX^-F)e^^ 
27r2ra4 

which means the exponential falloff of the free energy vs. the norm F of the strength tensor 
Fafj. On the other hand, we get the restriction for the upper limit of F, which follows from 
the requirement of the positiveness of the obtained free energy, F < 2g'^TXi. This means that 
the saddle-point configurations of the vector Sine-Gordon theory, we are studying, with too large 
strengthes are unstable. 

The other case is the degenerate one, when F — 0. In this case, 

7r2a4 ' 

which means that without the scale entering the system via the field x^, the free energy decays 
with the growth of loops as an exponent of the inverse square of their typical area. 
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